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INTRODUCTION 

It  is  well  known  that  the  strength  of  sea  ice  is  to 
a  large  extent  dependent  on  brine  volume,  and 
that  brine  volume  is  a  function  of  ice  salinity  and 
temperature.  Brine  volume  may  be  estimated  us¬ 
ing  the  empirical  expressions  of  Frankenstein  and 
Gamer  (1967)  or  the  theoretical  density-volume 
fraction  equations  of  Cox  and  Weeks  (1983).  The 
former  is  by  far  the  simplest  method,  because  the 
only  inputs  required  in  the  brine  volume  equa¬ 
tions  are  ice  salinity  and  temperature.  Since  the 
latter  approach  is  based  on  phase  relationships,  it 
allows  for  the  calculation  of  the  gas  volume  and 
thus  the  total  ice  porosity  and  the  brine-free  ice 
density.  With  this  method,  the  bulk  density  of  the 
sea  ice  is  also  needed  and  therefore  must  be  deter¬ 
mined.  In  principle  this  is  a  simple  task, 
but  it  can  be  difficult  to  do  in  the  field. 

Today,  ice  core  weight  is  measured  by  elec¬ 
tronic  scales  rather  than  triple  beam  bal¬ 
ances,  which  are  subject  to  uneven  support 
surface  and  wind  effect  errors  when  used 
at  unsheltered  sites.  Electronic  scale  weight 
measurements  can  be  highly  representative 
of  the  in-situ  ice  if  brine  drainage  has  not 
occurred  during  coring  or  after  the  ice  core 
was  removed  from  the  bore  hole.  Deter¬ 
mining  ice  core  volume  is  a  bit  more  trouble¬ 
some.  Ice  core  length  should  be  measured 
along  the  center  axis  using  a  special  caliper 
or  gauging  device  (Kovacs  1993).  A  meter 
rule  should  not  be  used  for  this  purpose  (it 
often  is),  because  the  ends  of  the  ice  sample 
may  not  be  parallel  to  each  other.  The  more 
difficult  measurement  is  ice  core  diameter. 

Again,  a  ruler  should  not  be  used.  A  pie 
tape  (Kovacs  1993)  or  a  caliper  with  long. 


wide  tongs  is  generally  satisfactory,  provided  the 
core  is  ridge-free  and  of  uniform  diameter.  Where 
this  is  not  the  case,  the  diameter  measured  will 
not  be  accurate  and  the  volume  determination 
will  be  off.  This  is  illustrated  in  Figure  1,  which 
gives  a  comparison  between  brine  volume  deter¬ 
minations  made  on  --10-cm-long  sections  from 
several  first-year  Arctic  sea-ice  cores  using  the 
equations  of  Cox  and  Weeks  and  Frankenstein 
and  Gamer.  The  regression  curve  through  the  data 
in  Figure  1  did  not  include  the  outlier  (+).  The 
outlier  indicates  where  a  measurement  error  oc¬ 
curred,  perhaps  because  of  a  volume  measure¬ 
ment  error. 

It  is  recommended  that  when  the  Cox  and 
Weeks  (1983)  equations  are  used,  and  they  must 
be  for  determining  gas  volume,  total  porosity,  and 


Figure  1.  Comparison  between  the  Cox  and  Weeks  and  the 
Frankenstein  and  Garner  methods  for  calculating  sea-ice  brine 
volume.  The  data  point  represented  by  +  is  considered  an 
outlier. 


Figure  2.  Sea-ice  bulk  density  and  conductivity  vs.  floe  thickness. 


brine-free  ice  density,  the  Frankenstein  and  Gar¬ 
ner  (1967)  equation  also  be  used  to  estimate  the 
ice  brine  volume.  When  the  two  brine-volume 
estimates  are  off  by  more  than  2%o,  a  measure¬ 
ment  or  data  logging  error  has  occurred.  If  the 
problem  cannot  be  resolved,  the  ice  property  val¬ 
ues  calculated  with  the  Cox  and  Weeks  equations 
should  not  be  used.  It  should  be  noted  that  the 
Frankenstein  and  Gamer  equations  for  determin¬ 
ing  brine  volume  are  not  valid  below  -23°C.  Only 
the  Cox  and  Weeks  equations  should  be  used 
below  this  temperature. 

The  physical  and  electromagnetic  properties  of 
21  winter  Beaufort  Sea  ice  cores,  obtained  by  the 
author  on  various  field  programs,  were  deter¬ 
mined  for  each  10-cm  increment  of  ice  core.  Physi¬ 
cal  property  calculations  were  made  using  the 
mathematical  expressions  of  Cox  and  Weeks 
(1983).  The  electromagnetic  properties  of  the  ice 
were  determined  using  the  procedures  of  Kovacs 
et  al.  (1987)  as  modified  by  Kovacs  and  Morey 
(1987).  Examples  of  the  calculations  are  given  in 
Tables  1  and  2.  Also  assembled  were  the  authors' 
data  from  an  additional  23  Beaufort  Sea  ice  cores, 
where  only  the  salinity  and  temperature  of  the 
winter  ice  were  determined  for  each  10-cm  incre¬ 
ment  of  core.  For  these  data,  only  the  brine  vol¬ 
ume  of  the  ice  could  be  calculated  using  the  equa¬ 
tions  of  Frankenstein  and  Gamer  (1967). 

From  the  above  data,  a  number  of  winter  ice 
floe  property  trends  can  be  seen.  Ice  floe  bulk 
density  vs.  thickness  is  shown  in  Figure  2.  The 
bulk  density  is  seen  to  decrease  with  ice  floe  thick¬ 
ness.  For  the  first-year  sea  ice,  the  decrease  is 
associated  with  brine  drainage  and  growth  rate 


processes,  which  reduce  the  vol¬ 
ume  fraction  of  the  heavier  brine 
entrained  within  the  ice  (Kovacs 
1996).  The  lower  multiyear  ice 
densities  are  the  result  of  this  ice 
containing  proportionally  less 
brine  and  more  gas,  especially  in 
the  freeboard  portion,  which  is  for 
all  intents  and  purposes  low-den¬ 
sity  fresh  ice.  While  density  is  the 
property  that  a  material's  strength 
is  often  related  to,  as  with  the 
strength  of  snow  (Kovacs  et  al. 
1969),  it  is  not  used  for  this  pur¬ 
pose  in  sea-ice  mechanics.  Tradi¬ 
tionally,  sea-ice  failure  strength, 
elastic  modulus,  and  such  have 
been  referenced  to  the  brine  vol¬ 
ume  of  the  ice.  This  brine  volume 
referencing  is  still  the  norm  for  sea-ice  tensile  or 
flexural  strength  test  presentations.  The  reason¬ 
ing  is  that  the  strength  of  a  material  is  a  function 
of  its  solidity.  Therefore,  the  more  brine  volume 
there  is,  the  less  solid  the  ice,  and  the  weaker  it 
will  be.  This  method,  of  course,  does  not  take  into 
account  the  added  porosity  associated  with  the 
small  gas  inclusions.  Since  brine  volume,  salinity, 
and  temperature  are  interrelated,  brine  volume 
needs  to  be  correlated  with  these  two  parameters, 
as  shown  in  Figure  3.  This  figure  clearly  illus¬ 
trates  that,  at  any  given  bulk  salinity  Sp,  the  bulk 
brine  volume  Vg  of  an  ice  sheet  is  dependent  on 
the  average  ice  sheet  temperature  While  this 
presentation  is  enlightening,  it  does  not  serve  the 
user  very  well.  More  useful  is  the  presentation  in 
Figure  4,  where  the  relationship  provided  allows 
for  the  determination  of  the  bulk  brine  volume 


Figure  3.  Sea-ice  bulk  brine  volume  vs.  bulk  salinity 
with  average  icefloe  temperature  as  a  parameter. 


2 


Table  1.  Physical  and  electromagnetic  properties  of  1.23-m-thick  sea  ice  at  100  MHz. 
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VB  =  41.64Sg 

r^  =  0,995 
SD  =  1.73 


Figure  4.  Sea-ice  bulk  brine  volume  vs.  average  absolute  ice 
temperature  and  bulk  salinity. 


vs.  ice-sheet  average  temperature  and  bulk  salin¬ 
ity.  It  should  be  realized  that  Figure  3  is  an  end 
vierv  of  the  3-D  data  presentation  in  Figure  4.  In 
Figure  3,  the  points  with  no  vertical  tails  are  within 
one  standard  deviation  (SD)  of  the  brine  volume 
surface  shown  in  Figure  4.  Points  with  tails  lie 
between  one  and  two  standard  deviations  of  the 
brine  volume  surface,  which  is  located  at  the  tip 
of  the  vertical  line. 


TENSILE  AND 
FLEXURAL  STRENGTH 


Similarly,  ice  floe  flexural  strength  Of  may 
be  estimated  using  the  empirical  equation 
(Timco  and  O'Brien  1994) 

Of  =1.76e-^-®®^.  (4) 

These  expressions  could  be  used  as  fol¬ 
lows.  Ice  thickness  is  determined  by  drill¬ 
hole  measurement.  Then  from  eq  2  the  bulk 
salinity  of  the  ice  is  estimated.  Ice  tempera¬ 
ture  is  determined  next  at  20-25  cm  below 
the  ice  surface.  A  linear  temperature  gradi¬ 
ent  to  the  ice/water  interface  (~1.7°C)  is 
assumed,  and  the  average  ice  sheet  tem¬ 
perature  is  calculated.  The  bulk  salinity  and 
average  temperature  values  are  then  used 
in  eq  1  to  determine  the  bulk  brine  volume. 
This  value  can  now  be  used  in  eq  3  or  4  to 
determine  Of  or  Of,  respectively.  It  should 
be  noted  that  Gavrilo  et  al.  (1995)  com¬ 
pared  various  published  expressions  for  deter¬ 
mining  the  flexural  strength  of  sea  ice  and  found 
that  eq  4  gave  0f  values  2  to  4  times  higher  than 
other  research  results.  Further  evaluation  of  the 
various  Of  equations  may  be  in  order. 

The  above  procedure  still  requires  a  drill-hole 
thickness  and  a  near-surface  ice  temperature  mea¬ 
surement.  These  measurements  may  not  be  nec¬ 
essary.  Using  the  formulations  for  determining 
the  electromagnetic  properties  of  sea  ice  (Kovacs 
and  Morey  1987),  the  low-frequency  (DC  to  -100 
MHz)  bulk  conductivity*  Og  of  a  model  sea-ice 


The  equation  for  predicting  the  bulk  brine  vol¬ 
ume  Vg  in  %o  as  a  fxmction  of  Sg  and  the  absolute 
value  of  (Fig.  4), 


Vb  =  41.64 


(1) 


along  with  the  expressions  of  Kovacs  (1996) 
for  estimating  the  bulk  salinity  (%»)  of  first- 
year  sea  ice  vs.  ice  floe  thickness  Tp  (cm) 

(Fig-  5), 


*  Sea-ice  conductivity  is  approximately  constant  between  DC 
and  100  MHz  (see  Tables  1  and  2). 


Sb=  4.606 +  91.603  /Tp, 


(2) 


Of  =0.816-0.069  V 


0.5 


(3) 


offer  the  opportunity  to  make  other  ice  prop¬ 
erty  assessments,  which  in  turn  can  be  used  to 
estimate  ice  sheet  strength. 

For  example,  the  horizontal  uniaxial  tensile 
strength  Of  of  an  ice  sheet  could  be  estimated 
using  the  empirical  equation  of  Dykins  (1970): 


Figure  5.  Combined  Arctic  and  Antarctic  ice  floe  bulk 
sdinity  vs.  thickness  (from  Kovacs  1996). 


4 


Vb  =  9.02  +217.09/T2  +  528.48cJb 
r2  =  0.996 
S  D  =  0.98 
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Figure  6.  3-D  presentation  of  sea-ice  brine  volume  vs.  DC  conductiv¬ 
ity  and  temperature  with  salinity  as  a  parameter  for  an  ice  density  of 
0.911  Mglm^  (a)  and  a  boresight  view  from  the  above  temperature 
window  (b). 
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increasing  ice  floe  thickness  as  re¬ 
vealed  in  Figure  2  affects  all  other  ice 
properties.  It  turns  out  that  this  ef¬ 
fect  is  favorable  to  the  remote  sound¬ 
ing  of  sea  ice.  For  example,  Kovacs 
and  Morey  (1991)  and  Kovacs  et  al. 
(1996)  have  shown  that  a  portable 
electromagnetic  induction  soiinding 
system  can  be  used  to  measure  re¬ 
motely  an  apparent  bulk  conductivity 
from  which  an  estimate  of  winter 
Beaufort  Sea  ice  thickness  can  be 
made  (Fig.  7).  When  the  Vg  vs.  at 
DC  values  listed  in  Tables  1  and  2 
and  from  20  other  ice  floes  were  com¬ 
pared,  it  was  found  that  the  calcu¬ 
lated  bulk  conductivity  of  Beaufort 
Sea  ice  was  directly  related  to  Vg,  as 
shown  in  Figure  8.  Therefore,  in  prin¬ 
ciple,  it  should  be  possible  to  esti¬ 
mate  both  undeformed  winter  Arctic 
ice  floe  thickness  and  bulk  brine  vol¬ 
ume  remotely  using  a  conductivity 
measurement  system.  Before  the  lat¬ 
ter  is  possible,  the  relation  between 
the  field-  measured  apparent  bulk 
conductivity,  which  is  highly  affect¬ 
ed  by  the  underlying  seawater,  and 
the  bulk  brine  volume  of  sea  ice  floes 
of  various  thicknesses  will  need  to 
be  determined.  Conductivity  mea¬ 
surement  at  a  frequency  well  above 
50  kHz  is  likely  to  be  necessary  to 
separate  the  sea-ice  conductivity  con- 


sheet  was  determined  along  with  its  brine 
volume  Vg  vs.  ice  temperature  and  salinity 
An  ice  sheet  bulk  density  of  0.911  Mg/m^ 
was  used  in  the  determinations.  The  calcu¬ 
lated  results  are  graphically  shown  in  Figure 
6.  There  is  a  decrease  in  brine  volume  with 
decreasing  ice  temperature.  In  addition  to 
ice  conductivity.  Figure  6  indicates  that  ei¬ 
ther  ice  salinity  or  temperature  is  required  to 
estimate  Vg. 

The  complex  interrelationship  between  the 
model  sea-ice  temperature,  salinity,  brine  vol¬ 
ume,  and  conductivity  is  well  revealed  in 
Figure  6.  The  results  suggest  that  determin¬ 
ing  Cg  of  natural  sea  ice  remotely  will  not,  of 
itself,  allow  Vg  to  be  estimated.  Not  revealed 
is  the  effect  of  density,  which  was  fixed  at 
0.911  Mg/m^  for  this  example.  The  decrease 
in  the  bulk  density  and  conductivity  with 


Drill-hole-Measured  Snow  and  Ice  Thickness  (m) 


Figure  7.  EM-31  instrument  conductivity  reading  at  -10  kHz 
vs.  ice  floe  snow  plus  ice  thickness  (from  Kovacs  et  al  1996). 
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Figure  8.  Measured  sea-ice  floe  bulk  brine  volume  vs,  the 
calculated  bulk  DC  conductivity  of  the  ice. 


Figure  9.  Ice  sheet  horizontal  tensile  strength  vs.  bulk  DC 
conductivity  and  brine  volume. 


Figure  10.  Ice  sheet  flexural  strength  vs.  bulk  DC  conductiv¬ 
ity  and  brine  volume.’ 


tribution  from  the  seawater  conductivity  con¬ 
tribution.  A  dual  high  and  low  frequency 
system  or  a  wideband  measurement  system 
could  be  used. 

In  the  interim,  an  off-the-shelf  electromag¬ 
netic  induction  system  operating  at  10  kHz 
can  be  used  to  estimate  ice  thickness  (Fig.  8), 
and  the  bulk  salinity  of  the  ice  is  estimated 
using  eq  2.  After  determining  the  average  ice 
sheet  temperature,  as  previously  described, 
eq  1  can  be  used  to  calculate  the  bulk  brine 
volume.  Ice  strength  may  then  be  calculated 
from  eq  3  or  4. 

The  relation  between  sea-ice  tensile  or  flex¬ 
ural  strength,  brine  volume,  and  the  calcu¬ 
lated  bulk  DC  conductivity  of  natural  sea  ice 
is  shown  in  Figures  9  and  10.  These  figures 
are  the  result  of  using  the  equation  in  Figure 
6  and  eq  3  and  4,  respectively.  In  the  future, 
using  similar  equations,  an  algorithm 
could  be  added  to  a  portable  electro¬ 
magnetic  induction  sounding  system 
that  would  provide  both  ice  thickness 
and  strength  from  the  apparent  sea-ice 
conductivity  measurements. 


COMPRESSIVE  STRENGTH 


Estimating  the  compressive  strength 
of  sea  ice  is  different.  While  the  uniaxial 
tensile  or  flexural  strength  of  sea  ice 
does  not  reveal  a  significant  strain  rate 
dependence,  at  the  strain  rates  of  inter¬ 
est  for  most  ice  engineering  problems, 
the  compressive  strength  does.  In  addi¬ 
tion,  in  reporting  sea-ice  compression 
test  results,  the  recent,  and  correct, 
method  is  to  relate  strength  to  total  po¬ 
rosity:  brine  volume  plus  gas  volume. 
Recently,  Timco  and  Frederking  (1990) 
developed  empirical  equations  for  pre¬ 
dicting  the  large-scale  uniaxial  uncon¬ 
fined  compressive  strength  of  Arctic  sea 
ice.  Their  equations  are  based  on  small- 
scale  strength  tests  and  the  various  con¬ 
trolling  factors,  such  as  ice  structure, 
brine  volume,  porosity,  loading  direc¬ 
tion  (horizontal-vertical),  and  strain 
rate.  Their  strength  equation  for  hori¬ 
zontally  loaded  columnar  sea  ice  is 


l-((t)B/270f® 
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Table  3.  Field-measured  and  calculated  test  data. 


Ice 

thick.^ 

cm 

Avg 

temp.^ 

Bulk  salinity 

%o 

Bulk  brine  volume 

%o 

Bulk  porosity 

%o 

Strain 

rate^ 

s~^ 

Horizontal  compressive  strength 

Measured^ 

MPa 

Calculated 

MPa 

156 

10.6 

5.50^ 

5.202 

38.19 

36.53^ 

52.87 

51.184 

7.90e-7 

0.81 

0.96^ 

0.38^ 

0.51^ 

0.548 

154 

10.5 

5.54 

5.20 

38.44 

36.26 

53.10 

51.39 

l.lOe-6 

0.49 

1.03 

0.46 

0.57 

0.60 

120 

12.9 

6.08 

5.37 

34.18 

32.94 

49.20 

47.97 

1.27e-6 

0.66 

1.07 

0.52 

0.62 

0.67 

129 

10.7 

5.94 

5.32 

38.50 

32.04 

53.16 

51.64 

9.20e-6 

1.00 

1.62 

1.12 

1.10 

1.17 

152 

10.5 

5.57 

5.21 

38.50 

36.82 

53.16 

51.94 

1.42^5 

1.30 

1.81 

1.30 

1.26 

1.35 

166 

4.2 

5.34 

5.16 

70.67 

67.46 

82.66 

79.20 

9.00e-5 

1.82 

2.15 

1.86 

1.79 

1.60 

147 

5.9 

5.65 

5.23 

56.83 

54.36 

69.96 

67.28 

9.10e-5 

1.86 

2.36 

2.04 

1.95 

1.88 

154 

4.3 

5.54 

5.20 

69.91 

66.86 

81.96 

78.65 

1.33e-4 

2.10 

2.33 

2.09 

2.03 

1.83 

146 

11.8 

5.66 

5.24 

35.72 

34.23 

50.61 

49.12 

1.56e^ 

2.86 

3.11 

2.70 

2.71 

3.00 

130 

5.1 

5.92 

5.32 

63.26 

60.85 

75.86 

73.18 

2.16e4 

1.88 

2.66 

2.50 

2.45 

2.28 

161 

7.3 

5.42 

5.18 

48.89 

46.74 

62.67 

60.37 

3.40e^ 

3.53 

3.41 

3.08 

3.08 

3.17 

129 

3.4 

5.94 

5.32 

83.08 

79.84 

93.96 

90.54 

7.50e4 

2.62 

2.86 

3.13 

3.24 

2.77 

a  =  EXXON  Test  Data 

1  =  Timco  &  Frederking  eq  6  or  7 

2  =  Eq2 

3  =  Eql 

4  =  Figure  13's  equation 


5  =  Timco  &  Frederking,  eq  5 

6  =  Timco  &  Frederking^  revised,  eq  8 

7  =  Figure  18's  equation 

8  =  Eq  10 


where  =  horizontal  uniaxial  unconfined  com¬ 
pressive  strength  (MPa) 
e  =  strain  rate,  s~^ 

(t)p  =  bulk  ice  porosity  (%o). 

Timco  and  Frederking  compared  values  de¬ 
termined  from  eq  5  with  in-situ,  full  ice  sheet 
thickness,  horizontal,  uniaxial  unconfined  com¬ 
pression  tests  performed  by  Exxon  on  Beaufort 
Sea  ice.  The  Exxon  large-scale  test  data  comprised 
ice  thickness,  average  ice  temperature,  loading 
strain  rate,  and  the  measured  strength  (Table  3). 
Since  the  bulk  porosity  of  the  ice  sheets  tested 
was  not  provided,  the  authors  had  to  make  these 
estimates.  To  do  this,  they  first  included  addi¬ 
tional  data  in  Figure  11  and,  following  the  same 
approach  as  Cox  and  Weeks  (1974),  ran  two  linear 
regression  curves  through  the  data.  Only  slightly 
different  from  the  Cox  and  Weeks  equations  (Fig. 
11),  the  new  equations  are: 

Sb  =  13.4- 0.174  Tp  (6) 

for  Tp  <  0.34  cm;  and 

Sb  =  8.0- 0.016  Tp  (7) 

for  Tp  >  0,34  cm. 


From  these  equations,  Timco  and  Frederking 
determined  Sb  for  each  ice  sheet  tested  (Table  3). 
Then  the  bulk  density  of  the  sea  ice  was  selected 
to  be  0.907  Mg/m^.  With  these  determinations 
and  the  average  ice  floe  temperature  given  in 
Table  3,  the  expressions  provided  by  Cox  and 
Weeks  (1983)  were  used  to  calculate  the  bulk  po¬ 
rosity  of  each  ice  sheet.  Note  that  these  values 
were  not  provided  by  Timco  and  Frederking  in 
their  paper.  The  missing  bulk  porosity  values  were 
recalculated  along  with  the  bulk  brine  volume  for 
each  ice  sheet.  The  values  are  listed  in  the  left 
column  of  the  brine  volume  and  porosity  data 
columns  in  Table  3.  The  strengths  measured  by 
Exxon  and  calculated  by  Timco  and  Frederking 
are  compared  in  Table  3  and  shown  graphically  in 
Figure  12.  The  regression  curve  in  Figure  13  shows 
that  the  calculated  values  are  generally  higher 
than  the  measured  ones.  This  offset  may  be  due  to 
small-  to  large-scale  scaling  ambiguities. 

To  bring  eq  5  into  better  agreement  with  the 
Exxon  test  results,  the  bulk  salinity  and  brine  vol¬ 
ume  of  each  ice  sheet  were  determined  using  eq  2 
and  1,  respectively.  Next,  the  Sb/  T^,  and  (|)b  data 
from  Tables  1  and  2  and  others  referenced  were 
plotted  as  shown  in  Figure  13a.  A  boresight  view 
of  the  data  as  seen  looking  into  the  salinity  "win¬ 
dow"  is  shown  in  Figure  13b.  This  view  clearly 
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Figure  11.  Sea-ice  bulk  salinity  vs. 
floe  thickness  (after  Cox  and  Weeks 


1974). 


Figure  12.  Timco  and  Frederking  esti¬ 
mated  full  ice  sheet  horizontal  uncon¬ 
fined  compressive  strength  vs.  Exxon's 
field-measured  values. 


reveals  the  interrelationship  of  average  ice  sheet 
temperature  and  bulk  salinity  with  the  bulk  po¬ 
rosity 

Using  the  Sg  values  determined  from  eq  2,  the 
average  ice  floe  temperature  listed  in  Table  3,  and 
the  equation  provided  in  Figure  13,  (t)B  was  esti¬ 
mated  for  each  of  the  Exxon  ice  sheets.  These  Vg, 
Sg,  and  (t)B  values  are  listed  in  Table  3  as  the  right 
column  under  each  respective  heading.  Using  the 
new  (|)b  values  and  the  related  strain  rates  in  Table 
3,  Oc  was  recalculated  using  eq  5.  These  so-called 
modified  values  were  then  compared  (Fig.  14) 
with  the  Oc  estimates  made  by  Timco  and 
Frederking  (Table  3).  As  can  be  seen  in  Figure  14, 
using  the  above  estimates  in  eq  5  gives 
values  that  are  comparable  with  those  of  Timco 


and  Frederking.  However,  the  modified  method 
for  determining  gives  values  that  are  still  offset 
from  the  Exxon  field  test  results  (Fig.  15).  To  alle¬ 
viate  this  offset,  eq  5  was  revised  to  read 


_  a7’c0.22 


ac=37e' 


l-((l)”-^V270) 


0.5 


■0.65.  (8) 


This  revision  brings  the  values  (Table  3)  into 
better  agreement  with  the  field  determinations, 
as  shown  in  Figure  16. 

If  the  DC  conductivity  of  the  ice  could  be  de¬ 
termined  remotely  along  with  ice  sheet  thickness, 
then  Sg  could  be  determined  from  eq  2  and  Vg 
estimated  from  an  expression  similar  to  that  given 
in  Figure  6  and  (t)g  from  Figure  17.  The  data  in 
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Figure  17  are  from  Tables  1  and  2  and  the  sources 
cited  above.  This  figure  indicates  that  nonde- 
formed  Beaufort  Sea  ice  floes  have  a  gas  volume 
between  about  1  and  2%. 

The  Gc  values  calculated  using  eq  8  in  relation 
to  the  controlling  parameters  £  and  <j)g  can  be  an¬ 


alyzed  using  a  two-parameter  power  function. 
The  relationship  is  a  further  development  of  the 
Norton  (1929)  secondary  creep  rate  equation  con¬ 
taining  a  power  function  of  stress: 

e  =  Ba”, 
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Figure  13,  Icefloe  bulkporosity  vs.  average  temperature 
and  bulk  salinity  (a)  and  boresight  view  from  the  above 
bulk  salinity  window  (b). 


where  B  is  an  empirical  parameter  that  takes  into 
account  ice  structure,  temperature,  activation  en¬ 
ergy,  and  other  factors,  a  is  the  applied  stress,  and 
the  empirical  exponent  n  is  generally  found  to  be 
about  3  for  ice.  Rearranging  for  a,  the  above  equa¬ 
tion  becomes 

o  =  Bi 

where  B;^  =  The  former  expression  was  used 
by  Glen  (1955,  1958)  to  evaluate  the  creep  behav¬ 
ior  of  solid,  not  porous,  freshwater  ice.  Since  the 
ice  porosity  has  a  profound  effect  on  a,  the  above 
equation  must  be  modified  to  include  a  porosity 
parameter  (j)  as  follows: 

a  =  B2e^^”  <!>'"/  (9) 

where  82  =  and  m  is  an  empirical  exponent. 

Equation  9  was  evaluated  by  substituting  the 
Gc  values  obtained  from  eq  8  for  G  in  eq  9  and 
using  the  related  controlling  parameters  and  £ 
given  in  Table  3.  The  result  is  shown  in  Figure  18. 
Equation  9  is  shown  to  statistically  fit  the  data 
extremely  well  with  an  value  of  0.993. 

As  more  data  become  available,  the  failure  sur¬ 
face  and  therefore  the  variables  in  eq  9  may  be 
better  determined.  Even  though  the  failure  sur¬ 
face  fits  the  data  with  a  high  value  of  0.993,  the 
data  suggest  that  this  surface  is  not  a  straight  line 
as  viewed  from  the  strain  rate  window  in  Figure 


Figure  14.  Modified  Timco-Frederking  ice 
floe  strengths  vs.  their  original  values. 
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CJ c(M)  =  0.61 5  +  0.862  a  c  (Ex) 
r2  =  0.921 
S  D  =  0.25 


I  I _ I _ I _ I - 1 - 1 - 1 - 
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ac(Ex).  Exxon’s  Horizontal  Unconfined  Compressive  Strength  (MPa) 


Figure  15.  Modified  Timco-Frederking  ice 
floe  strengths  vs.  Exxon  s  field-measured 
values. 


ac(Ex).  Exxon’s  Horizontal  Unconfined  Compressive  Strength  (MPa) 


Figure  16.  Icefloe  strengths  calculated 
using  eq  8  vs.  Exxon's  field-measured  G^ 
values. 


Figure  17.  Sea-icefloe  bulk  porosity  vs.  bulk 
brine  volume. 


Figure  18.  Equation-8-determined  ice  floe  strengths  Figure  19.  Equation-8~determined  icefloe  cr^  strengths  vs. 
vs.  bulk  porosity  and  strain  rate  (a)  and  the  boresight  bulk  porosity  and  strain  rate  (a)  and  the  boresight  view  of 
view  from  the  above  strain  rate  window  (b).  the  above  from  the  strain  rate  window  (b). 


18b.  The  three  data  points  at  a  strain  rate  of 
about  are  shown  to  be  offset  from  their 
predicted  failure  surface  position  at  the  tip  of 
their  vertical  tails,  suggesting  that  the  failure 
surface  curves  downward.  This  is  not  surpris¬ 
ing  if  B  contains  a  nonlinear  viscosity  compo¬ 
nent.  In  this  case,  while  stress  is  still  a  func¬ 
tion  of  strain  rate,  it  is  not  directly  proportional 
(Mellor  1986),  at  least  at  strain  rates  between 
about  10"^  and  10“^  as  indicated  by  the  stress 
vs.  strain  rate  curves  for  fresh  and  saline  ice 
given  in  Sanderson  (1988).  The  Exxon  data 
fall  in  the  ductile-to-brittle  strain  rate  transi¬ 
tion  region  where  the  creep  power-law  may 
not  be  valid,  that  is,  £  ^  because  n  is  no 
longer  a  constant. 

A  statistically  more  representative  failure 
surface  for  the  data  is  shown  in  Figure  19a. 
This  curved  surface  fits  the  data  with  a  cor¬ 
relation  coefficient  of  =  1.000.  A  boresight 
view  of  the  data  from  the  strain  rate  window 
in  Figure  19a  is  shown  in  Figure  19b.  Based  on 
the  data  presented  by  Sanderson  (1988),  be¬ 


low  a  strain  rate  of  about  10“^  the  failure  surface 
should  be  linear,  a -3.  However,  the  existing  data 
are  insufficient  for  developing  this  trend.  In  the 
interim,  for  strain  rates  between  10”^  and  10“^,  the 
equation  given  in  Figure  18  is  preferred  because  it 
does  not  contain  a  second  constant.  The  equation 
for  the  failure  surface  when  Vg  is  substituted  for 
(t)B  in  Figure  18  is  =150.75  £^-^^  In  this 

expression,  Vg  implicitly  accounts  for  the  bulk  gas 
volume  or  the  porosity.  In  using  either  ex¬ 
pression,  the  only  parameters  needed  are  the  load¬ 
ing  strain  rate  and  the  bulk  ice  porosity  or  brine 
volume.  The  latter  may  one  day  be  estimated  from 
remote  conductivity  measurements  along  with  ice 
thickness  using  a  hehcopter-bome  electromagnetic 
induction  sounding  system  (Kovacs  and  Holladay 
1990). 

A  comparison  of  the  data  derived  from  us¬ 
ing  eq  8  with  the  horizontal  uniaxial  unconfined 
compressive  strength  data  of  others  should  be  of 
interest.  The  data  compiled  by  Sanderson  (1988; 
see  his  Fig.  4.8)  are  shown  in  Figure  20  for  the 
strain  rate  range  of  our  data.  In  this  standard 
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Figure  20.  Icefloe  values  calculated  using  eq  8  from 
the  study  compared  with  the  data  provided  by 
Sanderson  (1988).  The  data  were  fitted  with  a  single¬ 
term  power  equation. 


Figure  21.  Icefloe  values  calculated  using  eq  8  from 
this  study  compared  with  the  data  provided  by  Sanderson 
(1988).  The  data  were  fitted  with  non-log-log  power 
equations. 


presentation,  no  corrections  have  been  made  to 
the  data  to  account  for  differences  in  test  tech¬ 
nique,  ice  porosity,  grain  size,  or  crystal  orienta¬ 
tion.  This  is  not  unusual  since  some  of  these  pa¬ 
rameters  were  not  measured  or  provided  in  the 
original  source  material,  as  is  the  case  for  the 
Exxon  data  used  in  this  report.  For  this  reason 
there  is  considerable  scatter  in  the  data.  The  up¬ 
per  regression  curve  shown  in  Figure  20  fits  the 
data  from  Sanderson;  the  lower  curve  represents 
the  data  generated  from  eq  8.  The  small  amount 
of  data  from  this  study  coupled  with  the  un¬ 
knowns  mentioned  above  does  not  allow  a  defini¬ 
tive  statement  to  be  made  on  the  difference  be¬ 
tween  the  two  lines.  Clearly,  the  data  from  this 
study  fit  well  within  the  scatter  of  the  data  from 
Sanderson.  Both  sets  of  data  are  statistically  bet¬ 
ter  represented  by  the  nonlinear  curves  shown  in 
Figure  21.  Here  again,  the  failure  surface  appears 
to  be  a  non-power-law  function  of  the  strain  rate. 

The  Timco  and  Frederking  (1990)  study,  which 
resulted  in  the  formulation  of  eq  6,  is  based  on 
small-scale  test  results.  How  applicable  this  equa¬ 
tion  is  to  full-scale  ice  force  assessments  can  be 
questioned  on  the  basis  of  small-  to  large-scale 
scaling  uncertainties.  To  avoid  this  problem,  the 
full-scale  Exxon  data  were  directly  analyzed 
using  the  values  listed  in  the  right-hand  bulk  po¬ 
rosity  column  in  Table  3.  Applying  eq  9  to  analyze 
the  and  e  data  of  Exxon  vs.  the  above-related 
(t)B  values  (Fig.  22a)  gives  the  following  expres¬ 
sion  for  the  failure  surface: 


Oc  =  1970  (|)^-^^ .  (10) 

In  the  boresight  view  of  the  data  from  the  strain 
rate  window  in  Figure  22b,  the  "outlier"  with 
the  long  offset  tail  at  a  strain  rate  of  7.9  x 
is  very  apparent.  However,  removing  this  data 
point  and  reanalyzing  the  remaining  data  gives 
a  new  equation  whose  calculated  values  do 
not  vary  by  more  than  0.04  MPa  from  those  cal¬ 
culated  by  eq  10. 

The  form  of  the  expression  for  the  curved  fail¬ 
ure  surface  shown  in  Figure  17  does  not  fit  the 
new  data  set  as  well  as  the  power  type  equation 
in  Figure  22.  This  may  be  seen  by  comparing  the 
statistics  listed  in  Figure  22  with  those  in  Figure 
23.  The  boresight  view  of  the  failure  plane  from 
the  strain  rate  window  (Fig.  23b)  also  indicates  no 
fit  improvement  for  the  outlier  over  the  failure 
surface  shown  in  Figure  22b. 

Substituting  the  new  Vg  data  for  the  <t)g  data 
and  reanalyzing  as  before  gives  the  results  shown 
in  Figure  24.  The  expression  for  the  failure  surface 
is  now 

Oc  =523.24  e  0-32  (11) 

As  with  the  bulk  porosity  data,  applying  a  curved 
failure  surface  through  the  vs.  bulk  brine  vol¬ 
ume  and  8  data  does  not  result  in  a  statistically 
better  curve  fit. 

It  is  interesting  that  the  vs.  8  trend,  as  ob¬ 
tained  with  eq  8  and  from  Sanderson's  data,  ap- 
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e.  Strain  Rate  (s”  ’ ) 


Figure  22.  Exxon  ice  floe  values  vs.  bulk  porosity  and 
strain  rate  (a)  and  the  boresight  view  of  the  above  from  the 
strain  rate  window  (b).  Data  were  fitted  with  a  two-term 
power  function,  i.e.,  eq  9. 


e,  Strain  Rate  (s~^) 

Figure  23.  Ice  floe  values  vs.  bulk  porosity  and  strain 
rate  (a)  and  the  boresight  view  of  the  above  from  the  strain 
rate  window  (b).  Data  were  fitted  with  a  non-log-log 
power  equation. 


r2  =0.938 
S  D  =  0.2^ 


X  § 


10-7  10-6  10-5  10-4 

e.  Strain  Rate  (s"'* ) 


Figure  24.  Ice  floe  values  vs.  bulk  brine  volume  and 
strain  rate  (a)  and  the  boresight  view  of  the  above  from  the 
strain  rate  window  (b). 
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Figure  25.  Boresight  view  from  the  bulk  porosity  win¬ 
dow  in  Figure  22a. 


Figure  26.  Boresight  view  from  the  bulk  brine  volume 
window  in  Figure  24a. 


pears  to  favor  a  curved  failure  surface  (Fig.  17 
and  19),  while  the  trend  resulting  from  the  analy¬ 
sis  of  the  Exxon  vs.  e  data  does  not  show  a 
strong  preference.  Additional  full  ice  sheet  tests 
should  help  to  further  define  the  failure  surface, 
especially  outside  the  area  of  the  current  data. 
Within  the  range  of  the  current  data,  10"^s“^  <  e  < 
ICr^s”^  and  ~25%o  <  Vb  <  80  %o,  eq  11  well  repre¬ 
sents  the  Oj.  strength  of  sea  ice  as  determined 
from  full-scale  in-situ  tests.  This  is  equally  true 
for  eq  10  and  its  counterpart  in  Figure  23a,  in 
which  (1)b  is  a  parameter.  Nevertheless,  eq  10  is 
preferred  because  of  its  simpler  form;  it  does  not 
contain  a  secondary  constant.  In  addition,  out¬ 
side  the  range  of  this  field  data,  the  latter  equa¬ 
tion  predicts  an  unrealistic  zero  ice  strength  be¬ 
tween  strain  rates  of  about  10“^  and  10“®s“^. 

The  relationship  of  bulk  porosity  and  strain 
rate  to  sea  ice  strength  is  clearly  revealed  in  Fig¬ 
ure  25.  The  effect  of  <|)b  on  is  shown  to  be  signi¬ 


ficantly  greater  at  the  higher  strain  rates.  As  would 
be  expected,  the  same  trend  is  maintained  when 
Vb  is  substituted  for  (()b  (Fig.  26). 

If  one  assumes  that  n  and  m  are  3  and  -1  re¬ 
spectively,  then  eq  10  can  be  simplified  to 

Oc  =2.7x103  e^/ 3  (t)-^  (12) 

Gj,  values  calculated  from  eqs  10  and  12  were 
compared  and  found  to  cross-correlate  with  an 
=  1.000  and  a  standard  deviation  of  0.012  over  the 
range  of  the  Exxon  full-scale  test  data. 

By  rearranging  eq  9,  one  obtains  an  equation 
for  the  secondary  creep  of  sea  ice: 

e  =  B3o”r,  (13) 

where  B3  =  Bj”  and  w  =  -mn. 

To  take  into  account  the  creep  process  of  ice 
with  time,  eq  13  can  be  modified  to  describe  pri¬ 
mary  creep  as  follows: 

e  =  B4a”<t)^f^,  (14) 

where  parameter  B4  will  be  (and  the  magnitude  of 
exponents  n  and  m  may  be)  different  from  those 
in  eq  13,  and  X  is  an  exponent. 

For  a  more  complete  description  of  the  creep 
process  of  ice  under  a  complex  stress-strain  state, 
the  reader  is  referred  to  the  recent  work  of  Fish 
(1993). 

ICE-STRUCTURE 
INTERACTION  FORCE 

As  is  well  recognized  by  engineers  engaged  in 
estimating  sea-ice  forces  on  offshore  structures, 
many  variables  enter  into  the  calculations.  In  the 
following,  an  estimate  is  made  of  the  horizontal 
force  that  could  develop  between  a  cold  first-year 
sea-ice  sheet  and  an  offshore  structure  in  the  south¬ 
ern  Beaufort  Sea.  In  this  scenario  only  the  struc¬ 
ture  width  (90  m)  and  the  degree-days  of  seawa¬ 
ter  freezing  since  a  stable  ice  cover  formed  (5000) 
are  known.  From  the  work  of  Anderson  (1961), 
Wen  et  al.  (1991),  and  R.  Lewellen  (pers.  comm.), 
the  relationship  between  Arctic  ice  floe  thickness 
and  the  cumulative  freezing  degree-days  since  a 
stable  ice  cover  formed  was  determined  (Fig.  27). 
This  curve  in  Figure  27  is  based  on  sea-ice  field 
measurements  with  a  natural  snow  cover.  From 
the  equation  for  the  curve  in  Figure  27,  the  ex¬ 
pected  ice  thickness  for  5000  freezing  degree-days 
is  176  cm,  and  from  eq  2  the  bulk  salinity  of  the  ice 
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Figure  27,  Sea-ice  thickness  vs.  seawater  freezing  degree-days  after 
a  stable  ice  cover  has  formed. 


sheet  is  estimated  to  be  5.15%o .  We  assume  the  ice 
sheet  has  an  average  temperature  of 
“10°C.  This  is  a  reasonable  value  for  cold,  late- 
winter  first-year  sea  ice  along  the  Alaska  Beaufort 
Sea  coast.  The  bulk  porosity  is  found  to  be  52.20  %o 
from  the  equation,  <|)g  =19.37  +  36.18 
given  in  Figure  13a. 

The  peak  failure  stress  in  sea  ice  occurs  at  a 
strain  rate  of  about  10“^s“^  for  small-scale  ice 
samples  with  an  aspect  ratio  less  than  2.  As  the 
sample  size  or  the  volume  of  ice  imder  load  in¬ 
creases,  the  peak  stress  occurs  at  a  lower  strain 
rate  (Sanderson  1988,  Gavrilo  1995).  At  full-scale 
ice  structure  interaction  when  the  aspect  ratio  is 
greater  than  50,  the  peak  stress  appears  to  occur 
at  an  effective  strain  rate  around  10”^  s“k  Since 
the  structure  width  d  is  very  much  greater  than 
Tp,  plane  stress  conditions  apply  and  the  appar¬ 
ent  ice  sheet  velocity  V  at  this  strain  rate  is 
(Sanderson  1988) 

y=  8e0.44^^  (15) 

or  1.4  m/h.  Using  52.2%o  for  and  10“^  s"^  for  e 
in  eq  5  gives  1.6  MPa  for  the  ice  sheet  peak  hori¬ 
zontal  uniaxial  unconfined  compressive  strength. 
For  plane  stress  conditions,  the  indentation  force 
F  between  a  vertical  structure  and  the  level  ice 
sheet  may  be  estimated  using  the  reference  stress 
expression  (Sanderson  1988) 

F  =  1.15  dTj^a^.  (16) 

For  d  =  90  m,  Tp  =  1.76  m  and  =  1.6  MPa  as 
given  above,  the  global  ice  force  against  the  struc¬ 


ture  is  estimated  to  be  291  MN  or  =  3.24 
MN/m.  Since  most  estimates  for  the  in¬ 
tegrated  limiting  pack  ice  force  lie  be¬ 
tween  10^  and  lO^N/m  (Parmerter  and 
Coon  1973,  Rothrock  1975,  Hibler  1980, 
Nevel  1983  and  Croasdale  1984),  the 
pack  ice  driving  force  would  seem  to  be 
the  design  limiting-force  load.  However, 
since  the  pack  ice  driving  force  can  be 
distributed  across  a  floe  that  is  much 
wider  than  the  structure,  the  concen¬ 
trated  force  at  the  structure  can  easily 
exceed  the  stress  required  to  fail  the  ice. 
Therefore,  for  this  model  and  no  other 
factors  being  considered,  the  design 
force  would  be  the  calculated  value. 

The  constants  given  in  eq  15  and  16 
are  best  estimates  and  may  need  to  be 
revised  (Sanderson  1988).  TTie  estimated 
force  of  291  MN  is  considered  an  upper  limit  for 
the  conditions  given.  Indeed,  for  a  similar  ice  sheet 
thickness  and  structure  width,  the  ice  force  devel¬ 
oped  against  the  structure  due  to  nonsimultaneous 
ice  crushing  failure  appears  to  be  less  than  half 
the  above  calculated  indentation  force  (Wright 
and  Timco  1994). 


SUMMARY 

A  constitutive  relationship,  in  the  form 
Oc  =  ^2  (|)g ,  was  developed  for  predicting  the 

imconfined  compressive  strength  of  first-year  sea- 
ice  floes  as  a  function  of  only  two  parameters,  the 
applied  strain  rate  and  the  ice  floe  bulk  porosity. 
TTiis  equation  was  evaluated  using  full  sea-ice 
sheet  thickness,  horizontal  unconfined  compres¬ 
sion  test  data,  so  no  uncertain  small-  to  large- 
scale  scaling  laws  were  involved.  The  magnitudes 
of  the  parameters  in  the  equation  were  found  to 
be  Oc  =  2.7  X 10“^  (|)"^.  Using  this  relationship, 

an  example  is  given  for  estimating  the  design 
sea-ice  crushing  force  against  a  90-m-wide  off¬ 
shore  structure. 

Relationships  are  also  provided  for  determin¬ 
ing  ice  sheet  tensile  and  flexural  strength.  The 
relation  between  sea-ice  bulk  conductivity  and 
strength  was  explored.  First-year  sea-ice  thick¬ 
ness  can  currently  be  measured  using  surface  or 
airborne  systems.  It  is  predicted  that  one  day  ice 
sheet  conductivity  will  also  be  remotely  measured 
and  that  this  will  lead  to  the  capability  of  re¬ 
motely  estimating  sea-ice  strength. 

Two  widely  used  methods  for  calculating  the 
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brine  volume  of  sea  ice  were  discussed.  Both  meth¬ 
ods  were  found  to  provide  similar  results,  pro¬ 
vided  due  diligence  is  exercised  when  making 
the  required  supporting  measurements. 
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